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i|Page 1 , (eq^, (U)): The right sides of the equations for coefficients 

L and N should he preceded hy minus signs. 


- - — OTjjincroj — nr-genersi^ De^conBrue'retr~on“aja ordinary two^ 

dimensioml basis • Annula r area reductions of the order of two to one^ 
or of higher orders, may occur, making it necessary to include the 
effect of radial variation in stream-filament height and curvature in 
the meridional plane in a blade -to -blade solution- It appears feasible 
to obtain a more con^jlete analysis of the impeller flow by a reitera- 
tive calculation procedure along two different types of relative flow 
surfaces. A description of the two flow surfaces, called Sp (a blade- 
to-blade type of suiTace) and Sg (a meridional type of surface) , and 
the equations for calculating the supersonic flow on both of sur- 

faces can be found in reference 1. 
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As a first approximation to the complete flow description, the 
stream-fllament-thl.ckness variation to he used in the calculation of 
the Si surface of revolution may he dete rmin ed hy an axlsymmetrlc 
solution in the meridional plane. A method, for such an eixisymmetric 
solution is described and worked, out in detail with examples in refer- 
ence 2. A solution of this type gives the flow variations along the 
hlade span; however, the hlELde-to-hlade solution must he approx ima ted, 
from the axlsymmetrlc stpersonic solution, hlade circulation, and con- 
tinuity. A total of 160 hours of hand computing is required for 
100 grid points in the meridional plane. A solution of the type of ref- 
erence 2 is necessaiy to give the stream-sheet thickness and. curvature 
for the hlade -to -hlade (S]_ surface) solution considered in this report. 
The equations for the solution of the surface of revolution are 

presented in reference 3- 

The purpose of the subject report is to continue the work of ref- 
erence 3 hy applying the method to a specific problem in order that the 
practicability of the method and the time involved in computing a solu- 
tion ntay he determined. The calculation procedure used in the applica- 
tion of the method is outlined in detail. 

Once this hlade-to-hlade solution is obtained, a finer approxima- 
tion to the flow may he obtained hy computing the flow on several 83 
surfaces and hy then repeating the computation on S]_ surfaces which 
are no longer surfaces of revolution as described in reference 1 . 

The tip stream, filament of a constant -tip -diameter supersonic 
impeller was chosen for Illustrating the numierical technique inasmuch 
as it was originally proposed to separate the stream-filament-thickness 
effect from the stream-curvature effect. For Sp surfaces of revolu- 
tion other than the tip, the calculation time required will he increased 
slightly over the l /2 hour per grid point estimated here, since, for 
this case, it is necessaary to include stream-filament slope terms. It 
is felt that the additional calculating time required will not he 
excessive. 

For completeness, the equations given in reference 3 have been 
derived from a different approach for a fluid with negligible entropy 
change and a constant-modlfled-relatlve stagnation enthalpy. The equa- 
tions are given in f ini te- difference form for the general surface of 
revolution, and the computing procedure is detailed for the case where 
the surface -of -revolution curvature is neglected. 
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SYMBOLS 

The following symbols are used In this report : 

A area normal to flow 

App critical area (at M = l) 

a local speed of sound 

h static enthalpy 

coefficients in characteristic equations 
K constant 

2, orthogonal coordinates on mean scnrface of revolution 

M local Mach number 

n tinit vector normal to surface of revolution 

P total pressure 

p static pressure 

r radial distance from axis of machine (fig. l) 

T thickness parameter 

V absolute velocity 

W relative velocity 

z distance measured along axis of machine (fig. l) 

P relative flow angle ^ tan“l 

X ratio of specific heats, 1.4 

8 flow angle correction term, equation (B2) 

e Prandtl angle correction term, equation (B1) 

T] first characteristic family 

X slope of characteristic curve on mean surface of revolution 
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11 Macli angle, sia”^ i/M 

V Prandtl-Meyer angle 

I second ctiaracteristic family 

p mass density 

a slope of mean surface of revolution in meridional plane, 

tan"^ W /tf 
r' z 

T normal thickness of stream filament of revolution 

ijr stream function 

CD angular velocity of "blade 

Subscripts : 

a stagnation conditions 

i inlet 

Z , 9 meridional and circumferential components 

r,z ccraponents in r- and z-directions, respectively 

t tip 

1,2,3 points of characteristic net 

o two-dimensional solution 

+,- characteristic slope of 1st and 2nd characteristics 

Superscripts : 

* dimensionless values, lengths are normalized by the tip radiiis 

and velocity by the stagnation speed of sound 

_ vector velocity 

DEVELOPMEIRT OF GEHMEIAL EQUATIONS 

The flow equations for supersonic flow on a general surface of 
revolution with variable stream-filament thickness t = t(Z) (see fig. l) 
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can "be developed from the equation of continuity and the absolute vor- 
tlcity normal to the surface of revolution. 

For steady flow the continuity equation hecomes : 

•a(TpW^r) S(TpW^ 

— — + -sr- 

If the flow is isentropic such that 



P 


1 

= Kh^"l 


the continuity equation may he written as : 



= 0 


( 2 ) 


Assuming that the absolute vorticity of the flow field is zero, then 

Vx V = Vx(W + OiXr)= VXW + 2o>=0 (3) 


The equation for the vorticity that involves the velocities and 'Wqj 

on the surfane of revolution is needed to coniblne with the continuity 
equation. The vorticity normal to the surface of revolution is a func- 
tion of these surface velocities . The vorticity normal to the surface of 
revolution in vector notation is 

n • (V X W + 2oj) = 0 


This relation, \daen expended in terms of and using the method 

of reference 4 (p. 210, eq. 9), is 



(4) 


For convenience it is desirable to express the continuity and vorticity 
equations in tenns of the relative velocity W and the flow angle p : 


W.J = ¥ cos p 


¥(p =» ¥ sin p 
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Differentiating equations (2) and (4) and substituting the relations 
for the velocity components, the following equations result where 


cot^p. Sw . tan P cot^p SW . ^ „ Sp Sp 

W ^ 


(a^r sin a , sin a , S In T 


^ P. |W_l^ + ^+tanp^ + 
W 5T W rSv ST r5^ 


1 h 

W cos p \ 


+ 2o3l sin a = 


inlying the method of characteristics to equations "(5) and (p) 
according to reference 5 (pp . 38-45) , two families of characteristic 
lines are obtained. The two typea of characteristic lines will be 
denoted by the characteristic parameters q and i, where t] varies 
along a | -constant- characteristic line and | varies along an Tj- 
constant- characteristic line. The characteristic equations become: 


W tan 


1 Sw Sp , CD r sin aj.SlnTj^ X_ sin a cos p j. 

an p ^ “ §1 ■ ‘"L a2 cos P-sin p) 


1 N cos P-L sin p ^ ^ ^ 
J cos P-sin p) “ 


where the slope of the tj- constant- characteristic line is 


\ = tan (p + p) 

aZ 

S imil arly, along a i -constant-characteristic line : 

1 SW Sp ^ CD^r sin a S In T sin a cos p 

+ _ tan p + + — + 

W tan p Bq St] L cos P-sin p) 


1 N cos P-L sin p ^ 0 

J (A^ cos P-sin p) Bq 


2935 



2935 


NACA rar 2992 


7 


and the i -constant-characteristic line slope is 


X = — =« tan (p - p) 


( 10 ) 


These characteristic lines are shown in figure 2. 

The coefficients J, L, and N on characteristic equations (?) and 
(9) are; 



A 


L = 




2 cos^ P 

03 r + 


sin a + a' 


2 S In T 


hi 


N = 1 - 




fW sin p 


+ 20b; 


sin a 


W cos p 




( 11 ) 


Characteristic equations (?) and (9) are the same as equations (l2a) 
and (l2h) of reference 3 for a flow with uniform I and negligible entropy 
change , where I is the modified- relative stagnation enthalpy defined by 


I = h + 


i . I 


A solution of the flow field is obtained by replacing the character- 
istic differential equations with a series of difference equations which 
are solved in a step-by-step procedure. Assume two characteristics 
emanating from points 1 and 2 and intersecting at point 3 as shown in 
figure 2 where the flow properties are known, one of them, q, with a 
slope equal to tein(p+p)i and the other, |, having a slope tan(p-p)g • 
The coordinates Z3, defining point 3 are determined from equa- 
tions (8) and (lO) . Thus 


riCqJj-qJ^) => (n constant) ^ 


^2^‘^3'^2^ “ ^-2^^3“^2^ constant) 


( 12 ) 


After point 3 is established from equations (12), W3 and P3 are deter- 
mined from the simultaneous solution of the characteristic equations. In 
finite difference form, equation (?) along ti and equation (9) along | 
are 
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iL Hi (“S-Wl) - I»1 + 

X- sin a cos g 1 R cos fi-L sin P~ | „ 

r (X_ cos 3 -sin ^ J (X_ cos P-aln PTjj^ ^ 


3 In T 
hi 


(t] constant) 


(13) 


H, («3-W2) + Os-Sa) - P" IT ° 

X^, sin a cos p 1 W cos P-L sin 3~ | / , _ i \ _ q 
r(X^ cos P- 8 in p) ^ J (X+ cos p-sin P)J 2 ^ 


S In, T 

ha 


(I constant) 


(14) 

For the case where a characteristic intersects a houndary rather than 
another characteiristic as shown in figure 2, only one characteristic 
eq.uation is needed., because either W 3 or will he known at the 

intersection point. In the case of an inverse problem either the veloc- 
ity or blade turning is specified, on the blade surface, ahd. in a direct 
problem the blade angle is known at the point of intersection. 

The conputational procedure to determine the properties and loca- 
tion of point 3 involves an Iteration process. The initial solution at 
point 3 is determined with the coefficients of •> (^3"^1,2) ^ 

(W 2 -VJ 2 ^^ 2 ), and (P3“Pi^2) "based on the values at points 1 and 2. In 

succeeding approxinKitions to the values at 3,. the coefficients based on 
the preceding values at point 3 eire determined and averaged with the 
coefficients at 1 and 2, so that after the first approximation to 3, all 
coefficients subscripted 1 and 2 are coefficients averaged between 1 and 
3 and between 2 and 3. The calculation procedure for the values at 3, is 
repeated until the solution at point 3 converges within the accuracy 
desired and the process is carried downstream step by step. 


CHARACTERISTIC EQUATIOHS FOR TIP STREAM FUAMEKT 

In this report the general characteristics method described has been 
simplified for application to the tip section of a constant-tlp- 
diameter supersonic conpressor. The numerical analysis takes into 
account a variable stream-filament thickness but Ignores the effect of 
stream-filament curvature in the meridional plane since this effect will 
be small at the specified tip section. In the case of a tip stream 
filament the following conditions hold: 


2935 



KACA TN 2992 


9 


Z =« z 

r = = constant 

T = t(z) 

CT = 0° 




J 


( 15 ) 


When these conditions are applied to characteristic equations (7) and 
(9), the chaxacteristic eqtmtions for the tip section become: 


1 

Sw 


cot p tan p 5 In T ^ Q 

(16) 

W tan 



tan p-cot p d| 

1 


ap 

cot p tan p S In T ^ Q 

f 1 7^ 

W tan 

p 


tan p-f-cot p dt] 



and from equations ( 8 ) aad (lO) their corresponding slope equations 
become: 


^t(<P3-^l) =« 23 - 22 ^) (t] constant) 

rt(<P 3 -cp^) =» ^_ 2 ( 2 ^ 3 “^ 2 ) constant) 


(18) 


If it is noted that the first term in equations (16) and (17) is 
the Prandtl-Meyer angle dv, the characteristic equations become the 
usual plane -flow characteristic equations (two-dimensional solution) 
plizs a correction term for stream-filament thickness as seen from the 
following equations: 


^ , cot p tan p d In T 

5^ ~ tan jx-cot 6 5^ ” 

Sv ^ cot p tan p ^ In T ^ Q 

3ri[ Bq ~ tan ji+cot p Sri 


The plane-flow characteristic equations are simply 


51 - 31 - ° 


Sv ^ ^p 


S3 


0 


(19) 

( 20 ) 


( 21 ) 

( 22 ) 


The analysis of the flow at the impeller tip was carried out using 
equations (19) and (20) together with the Ames tables (ref. 6 ). 
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niiOSTRATIVE EXAMPLE 

The characteristics metlaod out3JLned in the preceding sections for 
solving the supersonic flow equations has been applied at the tip of a 
constant -tip -diameter supersonic impeller to determine the blade-to- 
blade flow properties. The blade shape at the tip is shown in figure 3. 
An estimated boundary-layer -thickness correction has been added to the 
blade profile varying linearly along the axial blade length from zero at 
the leading edge to 10 percent of the pitch distance at the exit. The 
actual profile is shown by the dotted line in figure 3. The stream- 
filament -thickness variation from inlet to exit is assumed to be the 
result of a solution as described in reference 2. This variation at the 
tip section is shown in figure 4, where it is noted that the thickness 
is reduced at the exit to 40 percent of the value at the inlet. In the 
example, this is the axial-stream-fllament -thickness variation that is 
applied at the constant tip diameter (thickness remains constant tangen- 
tially) ; hence the stream-sheet curvature in the meridional plane was 
neglected. The quantities known at the leading edge were inlet angle 
Pi, -67.4°; inlet-flow deflection angle, 7°; suction-surface direction, 
-66.4°; inlet-Mach number M^, 1.6563; inlet-stream-filament thickness 
1^, 0.05; relative velocity W*, 1.3308; pitch distance, 0.2167; and the 
axial-chord length, 0.625. Dimensionless quantities are denoted with an 
asterisk. 

Special procedures not discussed previously were required in some 
phases of the numerical analysis because of the nature of the inlet con- 
ditions prescribed in the exEui 5 )le and because of the leading-edge shock 
wave interaction with the characteristic lines. These phases of the 
calculations and the determination of the characteristic blade-to-blade 
solution will be discussed in the following sections. 

Entrance region . - The solution of the flow equations by the method 
of characteristics assumes that' the boundary conditions are known on an 
initial curve, in this case an inlet line, along which a series of 
points, such as 1 and 2 in figure 2, Eire selected from which to initiate 
the characteristics solution. In establishing this line, it is neces- 
sary first to investigate the properties in the inlet region of those 
characteristic curves which are known as expansion waves. In accordance 
with the coordinate system selected for this example, these curves are 
characteristic lines of constant | where is tan (P-p) . When the 

absolute value of the wave angle, p-p, is less than 90° so that the 
expansion waves from the leading edge are contained within the blade 
passage as shown by the dashed line in figure 5, the inlet conditions as 
prescribed prevail up to the blade leading edge. For such a case the 
characteristic solution is started by selecting a series of grid points 
along AB where the flow is known and applying the iteration process 
described earlier at points 1 and 2 for the solution at 3. The proce- 
dure is possible since, in addition to the initial conditions along the 
characteristic line, the flow direction on the two blade surfaces is 
known. 
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However, with inlet conditions of the magnitude selected for this 
example, it was found that the absolute value of the wave angle, 3-p, 

“ at the leading edge was greater than 90° and the waves fell outside the 

blade passages. This type of wave system and the one described previ- 
ously are characterized by the magnitude of the axial component of 
velocity. If the eixial velocity component is supersonic, the waves will 
be contained within the blade passage, and if the axial component is 
subsonic, the waves fall upstream of the blade passage. For the latter 
condition where the flow upstream of each blade is influenced by the 
CD waves issuing, from the neighboring blade, the usual prescription of 

S unifonn-flow conditions at the entrance to the blade passage can no 

longer be made and it becomes necessary to study the flow conditions in 
this region further. 

In order to establish uniform-flow conditions at Infinity upstream 
of the blade row for steady-state flow, it is necessary that no disturb- 
ance exist at infinity. Any wave system that is present in front of 
the blades must originate from the blades themselves in such a manner 
that the expansion waves produced are cancelled by the formation, of eq.ual 
compression waves. Since only one family of uvaves, | -constant- 
characteristic family, extends upstream, a simple wave system will sat- 
isfy unif orm-flow conditions at infinity. 

For plane flow, the condition for a simple wave is merely v = v(^), 
P = P(l)^ which satisfies eq.uation (22) and equation (21) with the solu- 
tion V-p equal to a constant. This method cannot be used with equa- 
tions (l9) and (20) because V s= v(|), p = p(|) gives with equation (20) 
that 3 =» 0, which is not satisfied. The characteristic equa- 

tions (19) and (20) become somewhat simpler to study in terms of W sin P 
€ind a parameter 


1 

T-1 

T = Th W cos p 

When these variables- are introduced, equations (19) and (20) become 


1 

8(W sin p) 

X ^(1^ T) 

(23) 

W cos p -tan 

p 


1 

B(W sin p) 

a(ln T) 

(24) 

W cos P -tan 

p 



In this case the condition W sin p = 1(^)7 T =3 T(|) does not 
violate equation (24) but satisfies it identically; however, equa- 
tion (23) is not satisfied in general because l/w cos p -tan p is a 
function of Y, T, and T, and T is a function of | and t]. Differ- 
entiation of equation (23) with t](Y » Y(|), T = T(|)), gives the result 
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hi h ( 1 \hr 

\W cos p tan n / Brj 


0 


Thus there Is the solution Y equal a constant, and hence from equa- 
tion (23) T is also a constant. 


For the exar^ile considered in this report, these conditions were 
not used and the actual flow was approximated hy setting Y » Y(|) , 

T = T(^) and integrating equation (19). It is estimated that the maxi- 
mum eixor involved in V-p is less than 0.07° when the approximate 
value of the coefficient of d(ln t) is used in equation (19). Equa- 
tion (20) is satisfied exactly. The result of the integration is that 
V-p is seme function of thickness T, so that at the line AC (fig* 6) 
the quantity V-p is a constant (t is constant) which fulfills the 
condition for a single wave system upstream of the blades. 


An exact and detailed analysis of the flow at the entrance region 
ABC was not included in the original scope of this report. The results 
of the flow confutations to he presented are principally concerned with 
the flow downstream of characteristic AB; however, to determine the 
natijre of the waves in the entrance region, rou^ calculations approxi- 
mating the flow were made assuming the coefficients in the characteristic 
equations to he proportioned linearly along CB from the value at C 
to the value at B (fig. 6) . At the hlade leading edge, the inlet-flow 
direction was -67.4° and the suction-surface direction -66.4° so that 
there is a 1° expansion wave at the leading edge. A schematic sketch 
(fig. 6) shows the indications of the numerical results which are given 
in table I. The primed letters on figure 6 indicate points on the hlade 
surface along corresponding characteristics. These results show that 
compressions as well as expansion are propagated upstream from AC even 
though expansion waves are generated from surface CB because of veloc- 
ity changes along the characteristics. (These changes would he absent 
in plane flow.) The calculations further indicate little change in 
velocity and flow angle at the entrance to the blade passage so that 
assuming W and p constant along AC would have been a suitable 
assumption in this case. 

If a detailed analysis of the entrance region flow is not required, 
the flow properties may be determined by a method similar to the one 
described in the following section for obtaining the first inlet char- 
acteristic line AB contained within the passage. The fact has been 
utilized that W sin p equal to a constant along a |-constant charac- 
teristic gives a sinple wave system upstream of the blade row where T 
is considered. Since it has been assumed that W sin p is constant, 
equations (23) and (24) show that T is also a constant. 


2935 



NACA W 2992 


13 


Now with W sin 3 and T constant, the first expansion wave AB 
can he obtained by a simple graphical trial-and-error method. The rel- 
” ative velocity W* is fixed and a trial value for 3 selected at A. 

By the procedure shown in appendix A, line AB is solved piecewise, 
that is, from A to D, D to E, and so forth, establishing W* and 
P at each point. Qdie value of p determined at B must agree with 
the value of the blade direction at B- The calculation is repeated 
from A to B until the value selected for p at A gives in the 
w graphical solution the same p at B as given by the blade direction 

o) at that point-. In this manner the inlet characteristic line (expansion 

wave) of figure 7 was obtained. The length of the segments AD, DE, and 
so forth establishes the size of the subsequent characteristic network: 
since W* and p are determined at these points for beginning the 
characteristic solution. 

A Mach number of 1.6563 was prescribed at A in this example, 
fixing the relative velocity at 1.3308. An inlet angle of about -65° 
was desired; after the trial-and-error solution, an inlet angle of -67.4° 
at the leading edge was found to give the ‘.orrect inlet expansion wave. 

Characteristic solution blade-to-blade - The flow properties in 
the blade passage were determined using the "‘inite difference form of 
characteristic equations (19) and (20) as shi wn in appendix B. The 
resultant characteristic network of approxlma ely 100 grid points, where 
the flow properties were determined, is shown >n figure 7. 

The iterative procedure described earlier . ased upon averaged coef- 
ficients in the characteristic equations required from one to four iter- 
ations (usually two) to converge the relative velocity within one in the 
fourth decimal and the flow angle within one in the first decimal. The 
numerical setup used in the solution, together with a sample calcula- 
tion, is shown in table IV. 

Table II gives a conparison between the first and last iterations 
of the characteristic equations (l9) and (20) for four of the grid 
points indicated in figure 7. 

A calculation based upon continuity was carried out at the inlet 
and exit to determine what the nmgnitude of the relative velocity would 
be at the exit if the stream-filament-thickness reduction had not been 
taken into account as in a two-dimensional characteristic solution. It 
was found that the velocity W* at the exit would be approximately 
1.71 as compared with about 1.33 (see fig. 10) if the stream-filament - 
thickness effects were included. At each grid point the characteristic 
solution, compared with the plane-flow characteristics solution, indi- 
cated a reduction in relative velocity and absolute value of the flow 
angle when the stream-filament thickness was included in the character- 
istic equations. 
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Shock -wave . - The presence of the blade In a supersonic stream 
produces a shock -wave at the blade leading edge. In this example the 
shock wave is caused by the 7° deflection req.uired of the flow in turn- 
ing from the inlet angle of -67.4° to the pressure surface direction of 
-60.4° at the inlet. 

Across the shock iTave^ the variation in the properties themselves 
no longer can be considered infinitesimal and continuous as required in 
the application. of the characteristic equations; consequently, it is 
necessary to consider what happens when the shock wave crosses a char- 
acteristic line. Figure 8 illustrates a typical case where a shock wave 
from the leading edge intersects a characteristic line. The calculation 
of the flow properties for this situation is treated as follows : by 

ignoring the existence of the shock wave, an ordinary characteristic 
solution is carried out for point 3 (fig* 8) . The intersection point 
between the shock wave and the ^-characteristic line is then determined, 
art^ a linear interpolation is made between 2 and 3 for the flow proper- 
ties at c. Knowledge of the flow properties ahead of the shock and the 
fact that the -shock wave is a 7° compression wave permit the calculation 
of the flow behind the shock at point b from the general shock rela- 
tions; the flow can also be determined from reference 7, as was done in 
this e vamp Tff . The values at b are used to continue the next charac- 
teristic point solution. Point 3 is needed only to interpolate for the 
values at c and no longer appears in the characteristic network. 

It has been assumed in the calculations that the angle change (7°) 
is constant through the shock from the leading edge until it reaches 
the suction surface of the preceding blade. For a weak shock the flow 
deflection through the shock is independent of waves of the opposite 
family, so that a shock wave through a field of varying stream-filament 
thickness is also of constant deflection unless it is intersected by 
waves of the same family. As shown in figure 7, the leading-edge shock 
wave (heavy line) is not intersected by waves of the same family; there- 
fore, the calculations assuming a constant deflection are computationally 
consistent. However, if a finer grid had been used, the shock would 
have coalesced with waves of the same family and an adjustment of the 
flow deflection would have been necessitated. This situation occurred 
with the reflection of the shock wave from the suction surface where it 
is shown in figxire 7 to coalesce with expansion waves from the suction 
surface. The reflected shock wave was weakened from a 7° to a 2° 
deflection by the expansions from the suction surface. 

Streamlines. - A stream function ilr which satisfies the contin- 
uity equation (l) may be defined as follows : 


Si|r 


rrpW cos p 


(25) 
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« -TpW sin ^ (26) 

The stream ll nes shown in figure 9 were determined by integrating the 
stream function (eq.s. (25) and (26)) along a characteristic line extend- 
ing across the blade passage- Thus, 



BESUITS AMD DISCUSSION 

For a subsonic -axial- velocity component at the inlet to the blade, 
the family of waves traveling upstream must be defined by a simple wave 
system in order to establish uniform-flow conditions- at. infinity. In - 
this example the |-cons-tant-characteristic family extends iq)stream, and 
a simple wave system in this region where T is a constant is defined 
by V-p eq.ual to a constant. Also, if the curvature terms in the char- 
acteristic equations are neglected, W sin p equal to a cons-bant along a 
I -constant characteristic approx im ates a simple wave system in the 
entrance region. 

The characteristic solution shows a reduction in the absolute valu? 
of the flow angle and relative velocity in the blade passage from the 
•two-dimensional characteristic solution when a stream-filament -thickness 
reduction from blade entrance to exit is considered in the character- 
istic equations. This result would be expected in supersonic flow since 
a flow-area reduction produces a decrease in velocity. The flow prop- 
erties obtained from the last iteration using equations (l9) and (20) 
and the first iteration calculation check well, indicating that one 
iteration calculation may be sufficient. It is estimated that a solu- 
tion neglecting stream-filament curvature terms ob-tained with equa- 
tions (19) and (20), which is faclli-tated by using the Ames tables 
(ref. 6), would require l/2 hour of experienced computing time to con- 
verge a characteristic point. 

G?he streamlines determined are presented in figure 9. As would be 
expected for supersonic flow, the streamline spacing increases from 
pressure to suction surface. On the pressure surface are given the 
ma.-yiTniTm ijr-values divided by ijr -maximum at the inlet; these values are 
determined by integrating along |-cons-tant-characterlstic lines- These 
values show that continuity in the passage checked within ±2 percent- 
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Figure 10 stiows the velocity distrihution obtained on the blade 
surfaces. The abrupt decrease in velocity on the suction surface at 
z 0.25 is caused by the impingement of the leading edge shock wave on 
the suction surface at this point. The shock wave, noticeably decreased 
in strength, is reflected to the pressure surface at z « 0.505 where a 
lesser decrease in relative surface velocity occurs. The arithmetical 
mean velocity shown by the dashed line is approximately the same at inlet 
and exit, an indication that the area normal to the flow is the same at 
entrance and exit. 

The local static to relative upstream total pressure ratio on the 
blade surfaces is shown in figure 11. The sharp rise in pressure at 
z « 0.25 and z « 0.505 is caused by the shock wave. The average 
pressure shown by the dashed line gives an indication of the type of 
pressure distribution that may be expected from pressure measurements 
at the inpeller outer wall. 

Ho sinplified method of correcting for the stream-filament thick- 
ness could be found due to the interdependence of ^ and W in the 
characteristic eq.uations. The thickness correction as it is used here 
is applicable to thin sheets, since the thickness effect is actually 
propagated along characteristic lines in the meridional plane. The 
result is a translation downstream of the thickness effect from the 
axial station where the thickness actually existed. For thin sheets the 
translation is negligible. It may be desirable to reiterate the solu- 
tion with a new stream-filament-thickness distribution based upon the 
most recent velocity distribution. For the solutions suggested here to 
determine T* other fluid effects such as vorticlty and secondary flows 
will cause the actual flow to deviate from the hypothetical flow; hence, 
any gain by further refinements in t* may be masked by the other fluid 
properties which have been neglected in the solution. 


CONCLUDING EEMARKB 

The method of characteristics has been used to solve the Supersonic 
flow equations; this method gives the blade -to-b lade flow properties 
along an arbitrary stream surface of revolution where the meridional 
thickness associated with the surface of revolution is determined by an 
axisymmetric solution. 

The results Indicate an appreciable difference between the flow 
properties calculated with the two-dimensional characteristic equations, 
considering ordinary two-dimensional flow, and those determined with the 
characteristics method demonstrated here. For this particular example 
the effect of the stream-filament reduction between the blade -pas sage 
entrance and exit was to reduce the relative velocity and the absolute 
value of the flow angle. In the entrance region it is shown that 


2935 



MCA IN 2992 


17 


ex^panslon waves propagated along characteristics hecoine compression 
waves because of the stream-filament-thickness variation. This effect 
of the stream filament will result in an inlet-angle change upstream 
from that determined for two-dimensional flow. 

Expressing the characteristic equations in temis of the Prandtl- 
Meyer veirlahle v enables the Ames tables to be utilized in the com- 
putation procedure. The solution indicates that first iteration cal- 
culations gave flow properties of sufficient accuracy. Each grid point 
takes approximately l/2 hour to calculate, depending upon the number of 
iterations required. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Chio, April 17, 1953 
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APEENDIX A 


miET CHARACTERISTIC CALCULATION 

The following conditions were prescribed, at the leading edge point, 
5,0 (flg« 7), for the calculation of the inlet characteristic: 

“ 1*6563 
^5,0 =■ 1*3308 


^ 5,0 

^ 5,0 


0.05 


-67.4'- 


where 25^0 “ 0*00, “ l'*^-’ ^^5,0 0.2167. The value for ^5^0 

was found, to give the correct p at the characteristic -line inter- 
section, 0,0, on the suction surface after a series of solutions similar 
to that shown below. Along this inlet characteristic, 

T* » constant = (T*p*¥* cos P)5^0 “ 0.008568 

where p* = and. 

Pa 


(W* sin 6 ) 5^0 ~ constant = -1.2286 

With T*/l* aJii (W* sin p)^ as parameters, a straight-line eq.ua- 
tion can be developed, from which the inlet characteristic curve, and 
flow properties may be determined. Let 

T*/t* = d =■ p*W* cos p (Al) 

and upon substitution of 

p*W* = 0.578688 ^ 
into eq.uation (Al) , it becomes 

A* \ 

.578688 cos p (A2) 
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CO 

CD 

cn 


If 


e =* W* sin p 


the substitution of eq.uations (A2) and (A3) into 

1 ss cos^P + sin^P 


(A3) 


results in the following relation: 

1 


A 


,0.578688 


♦ \2 
cr I 


= i 


W 


♦2 


(A4) 


Equation (A4) is the equation for a straight line as shown in fig- 
ure 12 with an ordinate intercept equal to l/d^ and a negative slope 

of e^/d^. Another relation between r — and W 


*2 


\0. 578688 


can be 


detennined from the energy relation. This relation is 

_2 

1 


^0.578688 




(A5) 


The nature of equation (A5) is also shown in figure 12. The intersec- 
tion of the curve representing eqmtion (A5) with the straight line 
gives the value of the relative velocity. Straight lines in agreement 
With the inlet boundary conditions are calculated for a series of 
selected grid points. Once the velocity is determined from a plot sim- 
ilar to figure l2i, 'the flow single can be calculated from equation (A3) . 


The final inlet characteristic calculation is shewn in table III. 
Grid points were selected at approximately five equal intervals as 
shown in figure 7. 
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APPENDIX B 

CHARACTERISTIC CALCULATION IN BLADE PASSAGE 

Characteristic Equations (19) and (20) 

The variables u and p in chai*acteristic equations (19) and 
(20) can he expressed as a two-dimensional V and p plus a correction 
term- Let 


^o ^ 


P = Po + ^ 

Then equations (19) and (20) may he written as: 

d(hQ+c) - d(Po+5) + Gd(lnT*) =0 (ti constant) 

d(hQ+£) + d(pQ+8) - F d(ln T*) = 0 (| constant) 


(Bl) 

(B2) 

(B3) 

(B4) 


where F and G are the coefficients of the thickness term in equa- 
tions (l9) and (20). Since for the two-dimensional characteristics 

d^Q - dPQ = 0 (i) constant) 
dVQ + dpQ = 0 (| constant) 

equations (B3) and (B4) when integrated along the characteristic curve 
hecome : 


^3 - 63 + J G d(ln T^) = 0 (ti 

constant ) 

(B5) 

€3 + F d(ln T*) = 0(1 

constant ) 

(B6) 

For trapezoidal integration between characteristic points 1 and 3 and 
points 2 and 3 and for averaged coefficients F and G, the resulting 
correction terms are: 

F(ln T* - In Tp - G(LiT* 

- In Tp) 

(B7) 
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F(ln T* - la T*) + G(ln T* - In T?) 

63 5 5 (B8) 

Wlien the tvo -dimensional solution calculated in eq.uatlons (B9) to (B12) 
is utilized, €3 and 83 at point 4,3 can he calculated as shown in 
table IV. This procedure is representative of the type of solution 
required for characteristic equations (19) and (20). 


Plane-Flow Characteristic Equations (21) and (22) 
(•Dj^ - Vg) - (Pi " ^3) = 0 (n constant) 

^^2 “ ^3^ + i^2 ~ ^5^ ^ ^ constant) 
Solving for eund 

(V^ + Vg) + (Pg - Pi) 

u c= 

3 2 

(Pl + P2) " (^1 " ^2) 

^3 “ 2 

Calculations at grid point 4,3 are shown in table V. 


(B9) 

(BIO) 


(BU) 

(B12) 
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TABLE I. - APPROXIMATE EHTRANCE REGION FLOW PROPERTIES 


^■Point 

Relative 
flow angle, 

P, 

deg 

Relative 

velocity, 

¥* 

®'Point 

Relative 
flow angle, 

PS 

deg 

Relative 

velocity, 

¥*' 

C 

-66.4 

1.347 

C> 

-66.4 

1.; 

547 

S 

-67.0 

1.337 

S' 

-66.0 



R 

-67.3 

1.331 

R' 

-65.4 

t 

3 

Q 

-66.9 

1.337 

Q’ 

-64.0 



A 

-67.4 

1.331 

B 

-63.4 

1.; 

574 


^•Points defined on fig. 6. 


TABLE II. - CC»dPAEISON OF CHARACTERISTIC SOLDTIONS 
AT FOUR GRID POINTS 


Grid 

point 

Grid 

property 

Characteristic 
equations (l9) 
and (20), 
first iteration 

Characteristic 
equations (l9) 
and (20), 

n^^ iteration 


2,1 

P 

-59.82° 

-59.72° 



¥* 

1.4124 

1.4150 


3,1 

P 

-60.60° 

-60.54° 



¥* 

-1.4009 

1.4033 


4,3 

P 

-46.56° 

-46.61° 



¥ ♦ 

1.5539 

1.5527 


9,4 

P 

** * 

-11.66° 

-11.65° 



¥ * 

1.3864 

1. 3869 
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TABLE III. - INLET CHARACTERISTIC CALCULATION 



® 

© 

© 

® 

® 

® 

® 

Grid 

point 

z* 

arUitrax- 

ily 

selected 
along 
inlet 
line, AB 

T*, 

figure 4 
and @ 

a, 

0.008568 

© 

1 

d2' 

e2 

d2^ 

1.509 X 0 

1 

U*2’ 

asEvnne 
arbitrar- 
ily to 
find 
straight 
line 

Straight 
line 
points , 

0 - 0© 

4,0 

0.0085 

0.04875 

0.17575 

32.3749 

48.8699 

0.575 

.550 

4.2747 

5.4965 

3,0 

0.0155 

0.04777 

0.17936 

31.0849 

46.9227 

0.575 

.550 

4.1043 

5.2774 

2,0 

0.0220 

0.04690 

0.18269 

29.9620 

45.2276 

0.525 

.550 

6.2175 

5.0868 

1,0 

0.0270 

0.04618 

0.18553 

29.0516 

43.8534 

0.525 

.555 

6.0286 

4.9322 

0,0 

0.0318 

0.04555 

0.18810 

28.2632 

42,6633 

0.525 

.550 

5.8650 

4.7984 



® 

® 

0 

@ 

@ 

@ 

@ 


1 

W*, 

sin 3, 


M, 











Grid 

inter- 


-1.2286 

sin"^ 

ref. 6 

sin ^ ^ , 

@ - 

point 

section 


© 




deg 


fig. 12 


deg 


deg 


5,0 


1.3308 


-67.400 

1.6563 



4,0 

0.5571 

1.3398 

-0.91700 

-66.492 

1.6731 

36.705 

-103.197 

3,0 

0.5501 

1.3483 

-0.91122 

-65.675 

1.6901 

36.276 

-101.951 

2,0 

0.5432 

1.3568 

-0.90551 

-64.892 

1.7070 

35.861 

-100.753 

1,0 

0.5361 

1.3658 

-0.89955 

-64.099 

1.7251 

35.428 

-99.527 

0,0 

0.5297 

1.3740 

-0.89418 

-63.411 

1.7418 

35.038 

-98.449 
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TABIE IV. - CHAEACTERISnC CALCUIATION EQUATIOKS (19) AND (20) 




© 

© 

© 

© 

© 

© 

Grid 

point 

Grid 

point 

1 

P 

Pi 


cot (^2 

tan (2) 2 

©X(3) 

© - ® 

3 

2 

H 

^2 

cot (2) 2 

tan (Dg 

© + (D 

4,3 

4,2 

-56.00 

31.77 

-0.67450 

0. 61930 

-0.41772 

1.29280 


0,3 

-45.98 

27.04 

-.96636 

.51040 

-.49323 

-.45596 



© 

© 

© 

@ 


© 

Grid 

point 

3 

G, 

F, 

© 

© 

Table V 
colximn values 

®X(D 

©-© 

®+® 


cot [ 4 ] 

tan [a] 

Fo^ 

i 

4,3 

-0.32311 

1.08174 

-0.93546 

0.49988 

-0.46762 

1.43534 

-.43558 

-0.32579 

1.07356 



© 



@ 


@ 

Grid 

point 

3 

G, 

^3 

■'■3 


In @ 

In @ 

F, 

Average, 

®+@ 

^2 

2 

4,3 

-0.32445 

1.07765 

0.1368 

0.03430 

0.03877 

.03519 

-3.37260 

-3.25010 

-3.34698 




@ 



% 

(g) 

Grid 

point 

3 


@X@ 

S' 

@2-<g)i 

S' 

®2+@i 

Table V 
colunm values 

S' 

[3] 

S' 

B 

2 

X57.2958, 

deg 

2 

X57.2958, 

deg 

4,3 

-0.12250 

-.02562 

0.03974 

-.02761 

-1.93 

0.35 

32.66 

-46.91 



© 

0 

0 

0 

0 

0 

0 

0 

Grid 

^3' 

P 3 ' 

M, 

a 

^3' 

S' 

^+3' 

S3' 

point 

3 

0 

0 

(g) and 

s' 



0 

0 


+ 

+ 

ref. 6 

and 

X 

and 

+ 



0. 

©' 


ref. 6 

0 

ref. 6, 

0' 

0' 


deg 

deg 




deg 

deg 

deg 

4,3 

30.73 

-46.56 

2.1615 

0.7185 

1.5539 

27.02 

-19.54 

-73.58 
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TABLE V. - PLAHE-FLCW CHAEACTERISTIC CALCULATION 
EQUATIONS (21) AND (22) 


Grid Grid 
point points 
3 1 




^1 

Pi 

V2 

P2 

23.57 

-56.00 

31.73 

-45.98 


Pi |[(^i+©2 [©i+(D2 

+©2 -®i3 -®i+®2] 

•7-2 -r 2 


32.66 -46.91 


Grid 

point 

3 


© 

© 

a/Sa, 

^3" 

(^and 

©X© 

ref. 6 

0.7071 

1.5811 


ref. 6 




^-1-3" 

© + ©. 
deg 

0-0, 

deg 

-20.35 

-73.47 
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Figure 1. - General stream filament of revolution. 



Figure 2, - Blade section and cliaracteristlc 
calculating points on mean surface of 
revolution. 





Figure 5. - Expanalon waves 
contained In blade paaaage 
entrance section. 



Flgxire 6. - Wave system extending ^JpBt^eam 
in blade passage entrance section. 
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